Theory of electron spin decoherence by interacting nuclear spins in a quantum dot 
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We present a quantum solution to the electron spin decoherence by a nuclear pair-correlation 
method for the electron-nuclear spin dynamics under a strong magnetic field and a temperature high 
for the nuclear spins but low for the electron. The theory incorporates the hyperfine interaction, 
the intrinsic (both direct and indirect) nuclear interactions, and the extrinsic nuclear coupling 
mediated by the hyperfine interaction with the single electron in question. The last is shown to 
be important in free-induction decay (FID) of the single electron spin coherence. The spin echo 
eliminates the hyperfine-mediated decoherence but only reduces the decoherence by the intrinsic 
nuclear interactions. Thus, the decoherence times for single spin FID and ensemble spin echo are 
significantly different. The decoherence is explained in terms of quantum entanglement, which 
involves more than the spectral diffusion. 
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I. INTRODUCTION 

Irreversible processes of a microscopic system in con- 
tact with a macroscopic system are central to nanoscience 
and to quantum information science. A canonical exam- 
ple is the spin decoherence and relaxation of an elec- 
tron localized by an impurity, an electrical gate, or a 
quantum dot in a semiconductor, which have been ex- 
tensively studied both in theory^i^i'^i'^i'^i'^i'^i^i^i^°i"i^^ and 
in experimentsii^ii^iiSiiSiiLiSiiSiSfl Single dot spin relax- 
ation time Ti > 0.1 ms has been measured for differ- 
ent quantum dot systems at low temperaturefi24i^ii^ii& in 
good agreement with theoretical estimates AiS The sin- 
gle spin decoherence time T2 has a lower bound of 10 ns 
established by the measurements of the inhomogeneously 
broadened Tj* for either a spatial ensemble of many dots^^ 
or a time ensemble of a single dotii2ii2iS£ Spin dephas- 
ing by phonon-scattering in quantum dots is suppressed 
at temperature below a few Kelvins jiSiii leaving the nu- 
clear spins as the dominant mechanism for electron spin 
dephasing. Refs. 00 gave, respectively, a theory and a 
numerical study of the effects of off-diagonal electron nu- 
clear hyperfine interaction. Refs. .6j7,.21. gave treatments 
of the effects of nuclear dipolar interaction and calculated 
the electron spin decoherence time with an ensemble spin 
echo. 

Our pursuit of a quantum theory of decoherence with- 
out the restriction of a stochastic theoryi^ is motivated 
by the need to control the electron spin decoherence. The 
nuclear spins coupled to the electron by the hyperfine 
interaction are taken to be the sole source of decoher- 
ence as m Ref. m However, we treat the interaction 
between nuclear spins which will be shown to be impor- 
tant in the high magnetic field regime while the neglect 
of the nuclear spin interaction in Ref. 00 is valid in the 
low field regime. Our method of solution of the many- 
spin problem keys on the evolution of the two flip-flop 
states of each pair of nuclear spins and is thus termed 
the pseudospin method. We establish conditions for the 



validity of our method. The method is simple enough 
for many applications, including the more advanced de- 
sign of pulse control of the electron spin to eliminate the 
decoherence effects 1^ It also produces a simple physical 
picture, which greatly aids the applications. We divide 
the interaction between two nuclear spins into two types, 
intrinsic and extrinsic, respectively independent and de- 
pendent on the single electron spin state. The intrin- 
sic interaction consists of the dipole-dipole coupling and 
the indirect coupling mediated by virtual inter-band spin 
transitions via the hyperfine interactioniSiS^iSSiSSiSi The 
extrinsic nuclear interaction is mediated by virtual spin 
flips between each of the two nuclei and the single elec- 
tron due to the off-diagonal hyperflne coupling. 

The results presented here are (i) a basic solution of 
the decoherence dynamics, both for a single electron spin 
and for an ensemble of independent electrons; (ii) numer- 
ical evaluations for GaAs dots; and (iii) an analysis of 
the different time dependence of coherence for intrinsic 
and extrinsic nuclear-nuclear interaction. We will show 
that the extrinsic hyperfine-mediated nuclear interaction 
plays an important role in single spin FID. The spin 
echo not only refocuses the dephasing by inhomogeneous 
broadening in ensemble dynamics but also eliminates the 
decoherence by extrinsic hyperfine-mediated nuclear in- 
teraction. Thus, the decoherence times for single spin 
FID and spin echo are significantly different. The usual 
practice of inferring the single electron spin dephasing 
time from ensemble echo measurement could be prob- 
lematic. Note that in NMR literature, the difference in 
timescale of FID and echoes, when the inhomogeneous 
effect is excluded, was first recognizediS2i2S 

The electron spin decoherence arises out of quan- 
tum entanglement between the electron spin states |±) 
and the many-nuclear spin states \J'). When a co- 
herent electron spin state C+|-|-) -I- C_|— ) is prepared, 
the initial state of the whole electron-nuclear system is 
the product state, iC+\+) + C-\~)) (g) \J). In time t, 
the nuclear states associated with the two electron spin 
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states diverge, yielding an entangled state of the form, 
C+I+) \J+{t)) + C_|-) ® \J-{t)). The electron spin 
coherence is measured by \{J^ {t)\J^ {t)) \ when the nu- 
clear spin degrees of freedom are traced out. Our theory 
consists in a direct attack of the many nuclear spin dy- 
namics. The assumption of the pure electron spin deco- 
herence time T2 being much shorter than its longitudinal 
spin relaxation time Ti will be shown to lead to a simple 
effective Hamiltonian for the whole system of the form 
|±)_ff-|-(±| in terms of the nuclear spin Hamiltonians 

H±. This simplifies the electron spin coherence to the 
overlap of the two nuclear spin states, each following the 
evolution conditioned on one spin state of the electron, 
|(J|e-^-V^+*|^)|. 

This entanglement approach to decoherence has an in- 
teresting relation to the antecedents in the decoherence 
literature. The square of the decoherence is formally the 
same as the Loschmidt echo if H^, say, is regarded as 
with a perturbation, which is related to the decoher- 
ence of the nuclear spin system.SS A model of a single 
spin coupled to a transverse field Ising chain is used to 
study the effect of quantum phase transition on the deco- 
herence of the Ising chain.— A key change in the model 
could make it a study of the decoherence of the single 
spin in the Ising spin bath. 

We will show that the nuclear spin dynamics is dom- 
inated by the nuclear spin pair flips. The pairs can be 
treated as independent of one another and only the two- 
spin correlations need be taken into account in the in- 
teracting nuclear spin dynamics. On a time scale small 
compared with the inverse nuclear couplings but ample 
for the electron spin decoherence, the number of pair-flip 
excitations are small compared with the number of nu- 
clear spin pairs available for spin flip, which come from 
the randomization of the nuclear spin directions at a tem- 
perature higher than the nuclear spin temperature, i.e., 
10 mK ^ r < 1 K (still low enough to avoid the effects of 
the electron-phonon scattering). The cluster expansion 
by Witzel et al7 yields an equivalent pair-correlation ap- 
proximation. We will establish the pseudo-spin model in 
which the elementary excitations by independent pair- 
flips are just rotations of non-interacting 1/2 spins. 

Numerical evaluations for a GaAs dot then require no 
further approximation. The electron spin decoherence is 
calculated for a range of magnetic field strength (1-40 T) 
and various dot sizes. The pseudo-spin rotation yields 
the following analysis of the results. The e~* short-time 
behavior obeys n = 2 or 4 depending on the dominance 
of, respectively, the extrinsic and the intrinsic nuclear- 
nuclear interaction. In the long time limit, the crossover 
to the exponential decay [n — \) indicates the onset of 
the Markovian kinetics. 

The main body of the text gives a succinct account of 
the key points of the theory and the results of the compu- 
tation. Section defines the single electron system cou- 
pled to a bath of interacting nuclear spins. Section IIIII 
defines the electron spin coherence and formulates the 
quantum theory of its evolution. Section IIVI describes 



the pseudo-spin solution. Section IVI gives an evaluation 
of the decoherence for a GaAs quantum dot. Section IVTI 
serves as a brief summary of the main results. So as not 
to interrupt the flow of the essence of simple exposition 
of our decoherence theory, further details of the theory 
are grouped in the Appendices. 



II. SINGLE ELECTRON IN INTERACTING 
NUCLEAR SPIN BATH 

The system consists of an electron with spin vector 
So and N nuclear spins, J„, with Zeeman energies fJe 
and LJn under a magnetic field Bext, respectively, where 
n denotes both positions and isotope types (e.g. ^^As, 
^^Ga and ^^Ga in GaAs). The Hamiltonian of this sys- 
tem is described in Appendix ^ The interaction can 
be separated as "diagonal" terms which involve only the 
spin vector components along the field (z) direction and 
"off-diagonal" terms which involve spin fiips (see Ap- 
pendix 0. Because the electron Zeeman energy is much 
larger than the strength of the hyperfine interaction, the 
off-diagonal term is eliminated by a standard canonical 
transformation, with the second-order correction left as 
the hyperfine-mediated nuclear interaction. For the same 
reason, the off-diagonal part of the nuclear interaction 
contributes only when the terms conserve the Zeeman 
energies (so-called secular terms in the NMR terminol- 
ogy). Hence, the non-secular terms are negligible. The 
total reduced Hamiltonian is obtained from the transfor- 
mation in Appendix ^ for the limit of long longitudinal 
relaxation time (Ti 00), 



± 



\±)H+(±\ 



(1) 



with Hr. 



i^nJni ^'^d the interaction terms. 



H±=±Ha + Hb + Hd±He, (2) 

given by. 



, 4rje 

n<im 



(3b) 
(3c) 
(3d) 



where |±) are the eigenstates of S^, the summation 
with a prime runs over only the homo- nuclear pairs, 
the subscript A denotes the extrinsic hyperfine mediated 
interaction, B the off-diagonal part of the intrinsic in- 
teraction, D the diagonal part of the intrinsic interac- 
tion, and E the diagonal part of the contact electron- 
nuclear hyperfine interaction. The hyperfine energy;^ 
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determined by the electron wavefunction, has a typical 
energy scale En ~ a„ ~ 10^ s~^ for a dot with about 
10^ nuclei. The sum, A = J^n'^n, is the hyperfine con- 
stant depending only on the material. The intrinsic nu- 
clear spin-spin interaction has the near-neighbor coupling 
Bn.m ^ Dn,m ~ 10^ s^^. The hypcrfinc mediated inter- 
action, which is unrestricted in range and associated with 
opposite signs for opposite electron spin states, has an 
energy scale dependent on the field strength, An^m ~ 1- 
10 s-i for field ~ 40-1 T. 



III. DECOHERENCE THEORY 

The electron-nuclear spin system is assumed to be ini- 
tially prepared in a product state with the nuclear spins 
in a thermal state with temperature T, described by the 
density matrix 



P(0)-P°(0)®p^ 



(4) 



The time evolution of the reduced density matrix of the 
electron spin. 



p'{t) = TrNp(t), 



(5) 



obtained by tracing over the nuclear spins, may be ex- 
pressed in the form. 



(6) 



where p^, = (mIpI^^), and \p), H G {|+), |-)}. The 
superoperator or correlation function C^,v^^' can be ex- 
pressed in terms of the evolution operator and contains 
the information on the electron spin relaxation and de- 
coherence. 

The Hamiltonians of Eq. for the Ti oo limit 
conserves the electron 5*1 quantum number: [H^S^] = 0. 
Hence, the correlation function has following properties. 



^,At) - 1, (7b) 
£+^_(t) = C*_^+{t), (7c) 

and the specific expression for the free-induction decay. 



-ifl.t 



Tr^ 



(8) 



which can straightforwardly extended to dynamics under 
pulse control. 

The ensemble of nuclear spins, at temperature T > 
Bn,m,Dn^m,En, may be approximated by 
the density matrix. 



P 



(9) 



where \J) = 0^ j„ being the quantum number 
of nuclear spin n in the magnetic field direction. Pj is 
the thermal distribution factor. The correlation function 
C+,-{t) can then be generally expressed as, 

/:+,_(t) = ^Pje-^*-(*) \{J-{t)\J^{t))\ (10) 
J 

In FID, \ J^{t)) = e-*^**|J) and c^j{t) = {Qe + £j)t 
where £j = J^nJ^an is the contribution to the electron 
zeeman splitting from the Overhauser field in the nuclear 
configuration \J'). With a tt pulse to flip the electron spin 

at time r, we have |J^±(t > t)) = e-'^^^^'^'fe-'^^^lJ) 
and <l)j{t > t) = (f7e +£^)(2t - t). 

An important finding is that the coherence for the en- 
semble dynamics takes the factored form. 



C+.^{t) = C 



,_(t)x4"L(t), 



where 



/:;,_(i) = \{J-{t)\j+it))\ 



(11) 

(12a) 
(12b) 



J 



^+ - (t) characterize the electron coherence evolution in 
the single- system dynamics with the nuclear bath begins 
on a typical pure initial state IJ) = |j„). The single- 
system dynamics is solely determined by the spectrum 
of the elementary excitations from the initial state \J^), 
i.e., the nuclear spin pair- flips driven by Ha and Hb in 
Eqs. (|3al3b|) . For a sufhciently large number N of the 
nuclear spins, the spectrum is independent of the initial 
state with an error of the order of 1/^/N. Therefore the 
coherence factor £^ _ (t) can be pulled out of the sum 
over J' which results in the factorized form as in Eq. Hll|l . 
The independence of the nuclear pair excitation spectrum 
on the initial state has been numerically verified with 
a number of different initial states in a GaAs quantum 
dot with N = 10^. The independence of £^ _{t) on the 
initial nuclear state \J') makes the measurement of single- 
system coherence possible in principle, e.g., in single spin 
measurement with pre-determination or post-selection of 
local Overhauser field through projective measurementi^ 
The ensemble effect resides entirely in the factor 
>C?.^_(0, which may be read as the inhomogeneous broad- 
ening of the local field £j with distribution function 
Pj. The inhomogeneous factor (i.e., the distribution 
of the hyperfine energy £j) dominates the free induc- 
tion decay in the ensemble dynamics in the form of 
£f''_{t) = e-'"=*-(*/^2*)', with the dephasing time ~ 
VnA^^ ^ 10 ns as measurediA£*i2*i2iSli To single out 
the dynamical decoherence time from the Tj*, spin echo 
pulses can be applied to eliminate the effects of the static 
fluctuations of the local field. After a 7r-pulse applied at 
T, the inhomogeneous broadening part of the correlation 
function >C^''_(t) = 1 for i = 2r. The ensemble coher- 
ence peak at 2r is known as spin echo.^^ The spin echo 



4 



profile, i.e., the echo magnitude £+^_(2t) = £!^_(2r) 
plotted as a function of the echo delay time 2t, reveals 
the dynamical processes that leads to decoherence. 

IV. THE PSEUDO-SPIN SOLUTION 

The solution to the single-system evolution \J'^{t)) re- 
lies on the pair-correlation approximation explained here 
with more details in Appendix [21 and justified in Ap- 
pendix^ Within a time t much smaller than the inverse 
nuclear interaction strength, the total number of pair-flip 
excitations A'flip is much smaller than the number of nu- 
clei N . The probability of having pair-flips correlated is 
estimated in Appendix [Usl to be Pcorr ~ 1 - e"''^'''?/^ 
{q being the number of homo- nuclear nearest neighbors), 
which, as also shown by a posteriori numerical check, is 
well bounded by ~ 1% in the worst scenario studied in 
this paper. Thus, the pair-flips as elementary excitations 
from the initial state can be treated as independent of 
each other, with a relative error e < -Pcorr- Then the 
single-system dynamics \J'^{t)) can be described by the 
excitation of pair-correlations as non-interacting quasi- 
particles from the "vacuum" state driven by the 

"low-energy" effective Hamiltonian, 

^^-E^fe ^E^.^-'^feA (13) 

fe k 

which has been written in such a way that the pair- 
correlations are interpreted as 1/2-pseudo-spins, repre- 
sented by the Pauli matrix d^fc, with k labelling all possi- 
ble pair-flips, defined in more details in Appendix IC 41 
The time evolution from the initial state \J') can be 
viewed as the rotation of the pseudo-spins, initially all 
polarized along the -fz pseudo-axis: {^j, | tfe), under the 
effective pseudo-magnetic field, 

\i^ = {±2Au + 2Bk,Q.Du±Ek), (14) 

where, for the electron spin state |±), ±Ak and Bk are 
the pair-fiip transition amplitudes, defined in Eqs. |IC2|) . 
contributed by the hyperfine mediated coupling Ha and 
the intrinsic coupling Hb, respectively, and and ±£'fc 
are the energy cost of the pair-fiip contributed by the 
diagonal nuclear coupling Ho and the hyperfine inter- 
action He, respectively. The decoherence then can be 
analytically derived as 

^;,-(^) = Y{\{^k\^t)\ « n^"''^'' (15) 

k k 

where jV'^) are the conjugate pseudo-spin states. In FID, 

IV'^(i)) = e^*^fc*| tfc); while with a tt pulse to fiip the 

electron aXt^r, > r)) = e-'^?(*-^)e-*^^^| Tfc)- 

5\ = 1 — \ {'^k\'^k)\ possesses a simple geometrical in- 
terpretation: the squared distance between the two con- 
jugate pseudo-spins on the Bloch sphere, which quanti- 
fies the entanglement between the electron spin and the 
pseudo-spins. 
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FIG. 1: (Color online) Single-system FID for a dot with 
I/[ooi] = 2.8 nm and ro = 15 nm under various field strengths. 
The spin echo profile as a function of i = 2r for i?ext = 10 T 
is also plotted for comparison. The insets shows the field 
dependence of the FID decoherence time. 



V. RESULTS OF DECOHERENCE AND SPIN 
ECHO FOR A QUANTUM DOT 

In this section, we analyze the results of decoherence 
under free induction and spin echo conditions for a GaAs 
quantum dot. In numerical evaluations, the GaAs dot is 
assumed to have a hard-wall confinement in the growth 
direction [001] with thickness i^oQj^] and a parabolic con- 
finement with Fock-Darwin radius tq in the lateral direc- 
tions. The external magnetic field is applied along the 
[110] direction. For the indirect intrinsic nuclear interac- 
tion, we consider only the exchange partJ? The g factor 
of the electron^^ is taken as —0.13. The initial state \ J) 
is generated by randomly setting each nuclear spin ac- 
cording to a "high-temperature" Boltzmann distribution 
{Pj — constant). 

Fig. n shows the FID in single-system dynamics for a 
typical dot under various field strengths -Bcxt- The inset 
of Fig. n] shows the field dependence of decoherence time 
T2 which is defined as the time when the FID signal is 
1/e of its initial value. The strong field dependence of 
T2 demonstrate the significance of the extrinsic hyperfine 
mediated nuclear coupling until it is suppressed by a very 
strong field 20 T) . As shown in Fig. ^ the FID signals 
have significantly different decoherence times from the 
spin echo signals (See also Fig.^. In Fig.|3 we separate 
artificially the contributions from the extrinsic hyperfine- 
mediated and the intrinsic nuclear interactions and show 
their different dependence on time and dot-size. The spin 
echo profile is also plotted for comparison. 

A couple of justified simplifications can provide an un- 
derstanding of the effects of various mechanisms on the 
spin decoherence. First, the energy cost by the diago- 
nal nuclear coupling {Dk) can be neglected as it is by 
three orders of magnitude smaller than that by hyper- 
fine interaction {Ek). Second, for near-neighbor pair- 
fiips, the intrinsic nuclear interaction is much stronger 
than the hyperfine mediated one for the field strength 
under consideration. Third, for non-local pair-fiips, the 
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FIG. 2: (Color online) Separated contributions to the single- 
system FID by the extrinsic hyperfine-mediated nuclear cou- 
pling, the intrinsic nuclear interaction, and both, for a dot 
with L[ooi] = 6.2 nm and ro = 25 nm at Bcxt = 12 T. The 
spin echo profile with both mechanisms is shown in compar- 
ison. The inset shows the FID decoherence times resulting 
from only the hyperfine-mediated interaction (square symbol) 
or only the intrinsic interaction (triangle symbol) as functions 
of ro. 

intrinsic interaction is negligible due to its finite-range 
characteristic. Thus we can separate the flip-pairs into 
one subset, k S KA^ which contains 0(7V^) non-local 
flip-pairs, driven by the effective pseudo-magnetic fleld 

K. (±2ylfc, 0, ibiJfe) and a second subset, k G Kb-, 
which contains 0{N) near- neighbor flip-pairs, driven by 

« {2Bk,0,±Ek)- The conjugate pseudo-spins will 
precess along opposite directions in the non-local subset 
Ka, and symmetrically with respect to the y-z plane in 
the near-neighbor subset Kb- The decoherence can be 
readily grouped by the two different mechanisms as 

keKs keKA 

where hk = |h^| and sinc(a;) = sin(x)/x. We can see 
that the extrinsic hyperflne-mediated and the intrinsic 
couplings lead to the e"'*/-'"^ and the e~*^*/"^^'^^ be- 
havior in time shorter than the inverse pair-flip energy 
cost (which corresponds to the width of the excitation 
spectrum) , 

T2,B - h-^l''A-^I^N^I^\ T2,A « n,A-^N (17) 

where b is the typical value of near neighbor intrinsic 
nuclear coupling strength Bk and A = J2n the hy- 
perflne constant. Eq. (|17|l explains the dot-size depen- 
dence of the two mechanisms in inset of Fig. (21 When 
the two mechanisms are comparable, the single-system 

,2 

FID begins with e behavior and may cross over to- 
wards e~* decay as time increases, and this is actually 
observed in Fig. (21 As the time grows beyond the in- 
verse excitation spectrum width, the quantum kinetics 
becomes a stochastic Markovian process by building up 
the energy-conserving Fermi-Golden rule as indicated by 
the sine function in Eq. Hlt)|) . 
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FIG. 3: (Color online) Spin echo signal for dots of I/[ooi] = 
2.8 nm and various ro at _Bcxt = lOT. 

Fig. 01 shows the onset of the Markovian process by 
the crossover from the e~* short-time behavior to the 
long-time exponential decay e^*. 

In ensemble or repeated dynamicSfiiiSiSS FID will be 
dominated by the inhomogeneous broadening. In ESR 
experimentS(2^*^ electron spin flip pulses can be applied 
to eliminate the effects of the static fluctuations of the 
local fleld and spin echo proflles are measured. The echo 
decay time Th (deflned as the echo delay time at which 
the spin echo magnitude drops to 1/e of the zero delay 
value) is generally believed to give a quantitative mea- 
sure of the single-system FID time T2, as the direct mea- 
surement of the latter is of considerable difficulty with 
current experimental capability. We now show that the 
echo pulse will also modify the electron spin decoherence 
induced by the quantum pair-flip dynamics, and as a con- 
sequence, Th and T2 can be significantly different. As the 
electron spin is reversed by the 7r-pulse, the hyperfine- 
mediated transition amplitude Ak and the hyperfine en- 
ergy cost Ek for each pair-flip will change the sign after 
the pulse. Thus, the pseudo-spins driven by the extrinsic 
hyperflne-mediated nucear coupling (in subset Ka) will 
reverse their precession after the pulse and return to the 
origin at t = 2r, disentangling the electron spin and the 
pseudo-spins. So the decoherence driven by the extrinsic 
hyperflne-mediated coupling is largely eliminated in the 
spin-echo configuration (see Fig. |2J|. This phenomenon 
was first noted in the numerical simulation of a small 
system of ~ 10 nuclear spins in-^. For the pseudo-spin 
driven by the intrinsic coupling (subset Kb), the conju- 
gate pseudo-spins will switch their precession axis which 
also reverse the entanglement to some extent but no full 
recover can be obtained at the echo time. Finally, the 
electron spin coherence at the echo time can be derived 
as 

keKB 

Similar to the analysis for single system FID, the spin 
echo signal begins with e~(^'^/-^H ^ short-time behavior. 
The Fermi-Golden rule will build up as the time grows 
beyond the inverse pair-excitation spectrum width, which 
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FIG. 4: (Color online) Dot-size dependence of the decoherence 
times (see text), rH*" as circles, Th squares, T2 triangles, T2,b 
inverted triangles. is plotted (dashed blue line) for 

compare with T2,b- The dot thickness //[ooi] = 2.8 nm and 
Boxt = lOT. 



renders the quantum kinetics a Markovian process. The 
crossover towards the Markovian behavior becomes ob- 
servable when the width of the pair-excitation spectrum 
is greater than or comparable to the inverse of the initial 
dephasing time T^. For small quantum dots, the hy- 
perfine coupling and its gradient is larger, resulting in a 
broader excitation spectrum, and therefore the crossover 
behavior is observable, as shown in Fig. |51 For large 
quantum dots, where the pair-excitation spectrum is rel- 
atively narrower, the initial stage decay could already 
eliminate the spin coherence and thus the whole dephas- 
ing process could be described by the e~''^'^^'^H ) profile. 

The decoherence due to the intrinsic nuclear interac- 
tion is suppressed by the echo pulse as evidenced by 
the enhancement of the short-time decoherence time, 
T^h ^ y/2T2^B (see Fig. |21 and Fig. Ej). The onset of 
the Markovian crossover also manifests itself in the dif- 
ference between the initial echo time T^*^ and the overall 
echo decay time Th for small dots, as shown in Fig. ^ 



VI. SUMMARY 

In conclusion, we have presented the quantum the- 
ory for the electron spin decoherence by a bath of in- 
teracting nuclear spins under strong magnetic field and 
at finite temperature. Entanglement between electron 
and nuclear spins, established by their coupled evolu- 
tion, leads to the loss of electron spin coherence. The 
solution to the electron spin coherence amounts to solv- 
ing the many-body dynamics of the interacting nuclear 
spin bath, which are conditioned on different electron 
spin states. In the timescale of interest, the nuclear 
bath dynamics is dominated by pair-correlations among 
nuclear spins, which can be mapped into independent 
pseudo-spin excitations. Within the pair-correlation ap- 
proximation, the electron-nuclear spin dynamics reduces 
to the coupled evolution of the electron spin with the 
non-interacting pseudo-spin excitations, for which exact 



solutions are found. 

Decoherence behaviors in GaAs quantum dots are cal- 
culated as examples. We have demonstrated the signifi- 
cance of the extrinsic nuclear coupling mediated by the 
virtual electron spin flips, which manifests itself in the 
strong field dependence of the FID in single-system dy- 
namics. The calculated electron spin decoherence time in 
single-system FID varies from ~ 0.1 /xs to ~ 10 /is for field 
strength from 1 T to 20 T, and saturates as the hyperfine 
mediated coupling is suppressed by stronger field. The 
spin-echo pulse not only recovers the coherence lost by 
inhomogeneous broadening but also eliminates the deco- 
herence due to the hyperfine-mediated nuclear pair-flips 
and reduces the decoherence by the intrinsic nuclear in- 
teraction, leading to a spin echo decay time ~ 10 /is 
independent of field strength in ensemble dynamics. 

The theory presented here can be applied to other 
electron- nuclear spin systems, e.g., donor impurities in 
silicon, and may also be extended to even more general 
cases of a quantum object in contact with an interacting 
bath. 
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APPENDIX A: THE HAMILTONIAN 

The total Hamiltonian of the system of an electron and 
many nuclear spins is given by 

H = He + i?N + ^^cN + -f^NN, (A^l) 

composed of the single spin Zeeman energies Hq — ^^S^ 
and — X^n*-^""^" applied magnetic field along 

the z-axis, the hyperfine interaction iJoN and the intrinsic 
nuclear- nuclear interaction ifNN- 

The hyperfine interaction between the electron and the 
nuclear spins consists of the isotropic Fermi contact in- 
teraction and the anisotropic dipole-dipole interaction. 
The latter is negligible since the electron wave function 
in a quantum dot is dominated by the s-orbit states. The 
contact hyperfine interaction is given by^ 

HcN = ^a„Se-J„, (A2a) 

n 

an = ^7e7n^|^(R„)l', (A2b) 

47r 3 

where /iq is the vacuum magnetic permeability, R„ de- 
notes the coordinates of the 71th nucleus, 7„ is the nuclear 
gyromagnetic ratio, and 7e is the electron gyromagnetic 
ratio. It should be noted that while the effective (7-factor 
in the quantum dot determines the Zeeman energy Q,,, 
the free electron g-factor -1-2.0023 should be used for the 
hyperfine coupling.'^— 
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The intrinsic interaction between nuclear spins in- 
cludes the direct dipole-dipole interaction, the indirect 
interactions mediated by virtual excitation of electron- 
hole pairs, and the intra- nuclear quadrupole interaction. 
The direct dipole term is given bj^- 



semiconductor matrix of our interest and thus shall be 
the main ingredients of the intrinsic nuclear interactions. 

APPENDIX B: THE EXTRINSIC 
HYPERFINE-MEDIATED NUCLEAR SPIN-SPIN 
INTERACTION 



jjd 



E 



Mo Inlra 

47r i?3 



i?2 

n-.m 



(A3) 



with R„;m = R„ — Rm- The indirect nuclear interac- 
tion is mediated via the virtual excitation of electron-hole 
pairs by the hyperfine interaction between nuclei and va- 
lence electrons iSi^iSiSSiSi When the virtual excitation 
is caused by the Fermi-contact hyperfine interaction, the 
indirect interaction has the isotropic exchange form 



rrex 
-"NN 



(A4) 



named the pseudo-exchange interaction in the literature. 
The leading contribution of the pseudo-exchange inter- 
action for nearest neighbors in the host crystal can be 
expressed a.aS^^ 



ncx 



(A5) 



where is the effective gyromagnetic ratio determined 
by the renormalized charge density of the s-orbit elec- 
tron. The indirect exchange interaction has been exper- 
imentally studied by many researcher32^iSS421. 

When the virtual excitation of electron-hole pairs in- 
volves both the Fermi-contact and the dipolar hyperfine 
interactions, group theoretical analysis shows that the in- 
direct nuclear spin interaction has the dipolar form and 
is denoted as indirect pseudo-dipolar interaction in the 
literaturerS Lattice distortion can result in local elec- 
tric field gradients, inducing the intra- nuclear quadrupole 
interaction for nuclear spins with moment greater than 
1/2. The indirect pseudo-dipolar interaction and the 
quadrupole interaction are not included in the numeri- 
cal calculation in this paper due to the lack of experi- 
mental characterizations in the literature. Nonetheless, 
they can be well incorporated in our theory as contribu- 
tions to energy cost and transition amplitude for nuclear 
pair-flips [see Eqs. ljC2b|l and (jC2cp ] when reliable data 
is available. Furthermore, available studie o^^i^^i^^ show 
that the direct dipolar interaction together with indirect 
exchange interaction can mainly account for the broad- 
ening and lineshapes of NMR and NAR signals in the 



The contact hyperfine interaction can be separated into 
the longitudinal (or diagonal) part 



HcNj — anS^Jn, 



(Bl) 



and the transverse (or off-diagonal) part 

HcN.t 



EY(^eU:+^e-i+). (B2) 



The off-diagonal hyperfine interaction in the lowest order 
is eliminated by the standard canonical transformation. 



W = exp 



(B3) 



The residual second order term in the transformed Hamil- 
tonian H^cd = WHW~'^ is Ha in Eq. We neglect 

higher order terms whose effects are by at least a fac- 
tor of ^ (SleV^ /A)^"^ <C 1 smaller as compared to Ha- 
Note that while the zeeman energy D,,, ^ \Q — 100 ^eN 
for field - 1 - 10 T^ A/Vn - T^'^ - 0.1 /icV in GaAs 
fiuctuation dots.— 

The canonical transformation also rotates the basis 
states, |±) (g) \ J), \ J) = |j„), by a small amount 



W\±) ® \J) 



T ^wi±|T)®|j„±l) (g) |j;„),(B4) 

n ra^n 



where we have kept only the lowest order correction in 
terms of the expansion coefficient 



2 (fie - Un) 



\/j(j + l)-Jn(in±l), (B5) 



which is small number due to the inequalities, 17^ 3> 
S> a„. j = 3/2 for all three relevant isotopes ^^As, 
69Ga and ^^Ga. 



J 



The rotation of the state vector actually accounts for the rapid initial drop of the electron spin coherence. We 
show below that this reduction in the visibility of the spin beat is negligible under the field strength > 1 T. The state 
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evolution in the very initial stage {t < a^^^) can be expressed as 



-iHt 



|±) (g) \J) = W~^e~^"""'^W\±) ® \J) « iy-ie"*(^=+^"+-^<="'')*W^|±) \J) 



g±i(n,+£j-a;,.±a„)t\ I ±|2 



l±)»l^) 



_|_ ^ l^J^ _ g±i(Oe+f j-a;„±a„)t 



where £j — X^nin*^" the Overhauser energy of the 
electron spin under the nuclear configuration \J') which 
arises from the longitudinal (or diagonal) hyperfine in- 
teraction. To single out the process of visibility loss, we 
have omitted above the nuclear-nuclear coupling terms in 
iJjg(j^ The reduced density matrix of the electron spin 
can be derived as 

PlAt) « pU(0)[l-p+(t)]+P^_,_(OK(i), (B7a) 

P+{t) P-{ty 



with 



1 - 



(B7b) 



p±{t) =4:^{w^) sm . (B8) 



seed 




The interference between different frequency components 
on the RHS of Eq. ljB8(l will lead to a visibility loss of 
the coherence and the initial drop of the population with 



amplitude ^ p± ^ yflcvN/Aj ^ 1, occurring in 

the timescale ~ a^^. Therefore, in the high field limit 
(> 1 T), the electron spin flip by the nuclear spins is 
efficiently suppressed (Ti — )■ oo if phonon mechanisms 
excluded) . 

To conclude, the transverse (or off-diagonal) part of the 
hyperfine interaction has two effects: (i) the transforma- 
tion acting on the Hamiltonian results in an effective cou- 
pling between the nuclear spins {Ha in Eq. H3a|l ). which 
contribute to the pure dephasing of the electron spin co- 
herence; (ii) the transformation acting on the state vec- 
tor can be understood as a visibility loss. While the two 
processes coexist, we have shown that under the field 
strength > 1 T, the visibility loss is negligibly small. 
Therefore, the exact evolution of the electron nuclear sys- 
tem can be well approximated as 



-iHt 



±)<g>\J)^e~''''^^'\±)<E>\J) (B9) 



with the effect of the transverse hyperfine interaction well 
incorporated as the extrinsic nuclear coupling Ha- 



FIG. 5: (color online). The hierarchy structure for the nu- 
clear spin state evolution driven by pair-wise flip-lops. The 
k'l and k'^ pair-flips are possible only when the pair-fiips ki 
and ks have taken place, respectively. Some pair-flips occurs 
exclusively to each other, such as k2 and k^. 



flip-flop. Here, we examine the dynamics of the basic 
flip-flop process, based on which we build a hierarchical 
framework of many-spin basis states for dynamics, inves- 
tigate the spin correlations, and construct a pseudo-spin 
method. 



1. The basic nuclear spin excitation 

The transition for the pair-flip driven by the operator 
J+J~ consists in a pair of nuclear spins in state \jn,jm) 
goes to \ jn + l,jm — 1) if permitted, i.e., jn -I- 1 < j and 
jm — 1 > ~J for spin j. The shorthand k is used to de- 
note this pair transition. Such a fc-th pair-flip transition 
between two many-nuclear spin states is described as 



APPENDIX C: SOLUTION OF THE DYNAMICS 
OF THE ELECTRON-NUCLEAR SPIN SYSTEM 

The elementary process driven by the interaction 
Hamiltonian, Eq. is the pair- wise homo- nuclear spin 



\±)\J)^\±)\J,k). (CI) 

The transition matrix element ±Ak + Bk and the energy 
cost Dk ± Ek, the sign conditioned on the electron spin 
state |±), can be derived from Eqs.©, 
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Ak 

Bk 
Dk 



{J,k\HA\J) 



VJ (i + 1) - in Un + 1)VJ U + 1) - jm (jm ~ 1), 



{J, k\HB\J)^ S„,„VJ (i + 1) - .?« (Jn + 1) U + 1) - (jm " 1) 

A:| i^D 1^, k) - {J\ Hd \J) = Dn^n'in' - Dm,rn'jm' - Dn,^, 

n' rn' 

{J, k\ He \J, k) - {J\ He \J) - (a„ - a™) /2. 



(C2a) 

(C2b) 
(C2c) 

(C2d) 



2. The hierarchy of nuclear pair-flip states 



The evolution of the nuclear spin state \J^{t)) = 
g-zff ^\J) can be formally described by pair-flip transi- 
tions in a hierarchy of basis states. The hierarchy is com- 
posed of the seed state, \J), the first-generation states, 
each \J, k) generated from the seed state by the fcth pair- 
flip, the second-generation states, each \ J,ki,k2) gener- 
ated from the first-generation state \ J, ki) or \J, ^2) by 
the fcith or the fc2th pair-flip, respectively, and so on. A 
many-nuclear spin state can be expanded in this basis as 



fci 



in which the wavefunction satisfies the equation 



(C3) 



dtC 



JM 



-iE 



.k^....jkp Jj^ .k\....^}ip 



j=l kj 

tY iBl±Al)C^,^^ 



(C4) 



where the energy Ej f.^ ^.^ is the eigenenergy of the ba- 
sis state fci, . . . , kp) under the Hamiltonian H]j±He, 
respectively. The hierarchy description of the nuclear 
spin dynamics is illustrated in Fig. |S| 

The many-body nature of the problem lies in the fact 
that the pair-wise flip-fiops are correlated in general. The 
correlation between two pair-flips can be developed in the 
following three cases: 

1. Exclusive pair-flips as shown in Fig.|3(a) and (b). 
When one pair-flip (ki) has taken place, the other 
one (fc2) which shares one or two spins with the 
flipped pair cannot occur any more, and vice versa. 

2. Subsequent pair-flips as shown in Fig.|Hl(c) and (d). 
When two pair-flips share one or two nuclear spins, 
one pair-flip is possible only after the other one has 
already taken place. 



Ill II IIII 

(a) (b) (e) 



III li Il-il 



(c) 



(d) 



(f) 



FIG. 6: (Color online). Possible correlations between two nu- 
clear spin pair-flips (fci and k2, indicated by solid and dotted 
arrows, respectively), (a) and (b), the two pair-flips sharing 
one or two nuclei occur exclusively to each other, (c) and (d), 
one pair-flip can take place only after the other pair has been 
flipped, (e) The energy cost of one pair-flip depends on 
whether or not the other pair has been flipped, when they 
involve spins in neighborhood, (f) The two pair-flips are in- 
dependent of each other when they are far apart from each 
other. 



Neighboring pairs as shown in Fig.|3(e). Two pair- 
flips can be correlated even when they don't share 
a nuclear spin but are in neighborhood, since one 
pair-flip will change the nuclear spin configuration 
in the neighborhood of the other pair and thus mod- 
ify the energy cost {Dk) due to the diagonal nuclear 
spin interactions of one pair-fiip. This correlation 
can be clearly seen from Eq. (|C2cl) . 



When two pair-wise flip-flops are not in neighborhood 
[Fig. El (f)], they are independent of each other. 

With increasing numbers of pair-flips, the nuclear spin 
state involves basis states further up in the hierarchy 
structure depicted in Fig. [51 and higher order correlations 
may occur, making the solution of nuclear spin dynamics 
a formidable task in general. 



3. Pair-correlation approximation (PCA) 

Two pair-flips are correlated only when they are lo- 
cated in neighborhood. We can estimate the probability 
of having pair-flips correlated in the pth generation of 
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the hierarchy: ii p — I pair-flip excitations have been 
generated, the probabiUty of having the pth pair out- 
side the neighborhood of all the previous ones is about 
{l — q{p— 1)/N). By induction, the probability of having 
p pair-flips uncorrelated is 

1 - Pcorr ^ n (1 - - exp i-qpyN) . (C5) 

The number of pair-flip excitations at time t may be 
estimated by 



p=l 



(C6) 



If qNl^^{t) ^ in the timescale of interest, the probabil- 
ity, PcoiT, of having pair-flip excitations correlated is neg- 
ligibly small. Thus, by removing the few states reached 
via subsequently correlated pair-flips and adding few 
states containing exclusive pair-flips, the exact Hilbert 
space can be mapped into the tensor product of two- 
dimensional Hilbert subspaces, each of them correspond- 
ing to a pair- flip available from the seed state | namely 



{\J)AJ, h),\j,ki,k2),...} 



|{|U>,|U>},(C7) 



where the index k runs over all possible pair-flips from 
the seed state to the first generation. The mapping can 
be explicitly expressed as 



\J) 
\J, fci) 
\J, ki,k2) 



k 

ife.) (g) I 
k^ki 

iki) \ ik2) 



(C8a) 
(C8b) 
Tfe), (C8c) 



k^ki ,k2 



and the nuclear bath state can be factorized into ele- 
mentary excitations as 

=e-^J*(g)(.9±W| h)+fHt)\ U» • (C9) 



Furthermore, the energy cost of a pair-flip is assumed 
independent of whether or not another pair-flip has oc- 
curred in the neighborhood. So all the elementary exci- 
tations are treated independent of each other and their 
dynamics is determined solely by their own energy costs 
and transition matrix elements as 



(ClOa) 

{Bk±Ak)gt, (ClOb) 



with initial conditions 5^(0) = 1 and /^(O) = 0. 
The above approximation strategy is denoted as the 



pair-correlation approximation (PCA) where all pair- 
correlations among the nuclear spins are kept and higher 
order nuclear correlations are neglected. 

In PCA, we will calculate physical properties by using 
Eq. (jC9|) as the bath state at time t instead of the exact 
state Eq. I|('4p . In this paper, the electron spin coher- 
ence of interest is essentially the state overlap between 
two differently driven bath state as shown in Eq. (|12a|l . 
The relative amount of change to the bath Hilbert space 
structure by PCA is Pcorr- Therefore, the relative error 
of PCA in calculating the coherence is bounded by 



<Pc. 



1 



exp I 



(Cll) 



4. The pseudo-spin model 



The uncorrelated pair-flip dynamics given in Eq. ljC10|l 
is nothing but the independent evolutions of two-level 
systems {| tfc)^ ifc)}- Thus the electron- nuclear spin 
dynamics \ J^{t)) from the initial state \ J) is mapped, 
by PCA, to the precession of N non-interacting pseudo- 
spins under pseudo-flelds conditioned on the electron spin 
state 



h± = {±2Ak + 2Bk,Q,Dk±Ek), 



(C12a) 
(C12b) 



where = (o"^, ct^, ct^) is the Pauli matrix for the 
pseudo-spin corresponding to the fcth pair-flip available 
from the seed state \J). Note that by the Hamiltonian 
mapping from Eq. (0) to Eq. (|13|l , the energy of the seed 
state has been shifted by a constant [Dk ± Ek) /2, 
which does not affect the calculation of \ {J~ {t)\J^ {t))\. 
The directions (indicated by the superscripts x, y, and 
z) for the pseudo-spins are not defined in the real-space 
but in a fictitious pseudo-space. 



APPENDIX D: RANGE OF VALIDITY OF THE 
PAIR-CORRELATION APPROXIMATION 

Within the PCA, the averaged pair-flip number iVflip 
defined in Eq. (|C6|1 has a simpler form 



iVfi 



ip 




(Dl) 



with which we can make an a posterior self-consistency 
check of the validity of the approximation. If Nj\ip ob- 
tained from Eq. ljDl|) is small, we conclude that it also 
faithfully refiects the number of pair-flip excitations in 
the exact dynamics and therefore, the error estimation 
would be faithful and the physical property can be cal- 
culated based on PCA with a relative error bounded by 

Eq. icriTi . 

The total number of pair-flip excitations can be divide 
into two parts: iVflip {t) = A^'flip^A {t) + A'^Hip^ {t), where 
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■^flip.A is the number of non-local pair-flip excitations and 
-^flip.B is the number of local pair-flip excitations that 
have been created. A'fiip^A {t) and A^fiip^B (t) have very 
different dependence on time and system parameters, and 
we analyze them separately. 

In free-induction evolution, the number of non-local 
pair-flip excitations is given by, 

where Ma ^ N'^ is the number of non-local nuclear 
spin pairs. Since the evolution of the non-local pair- 
correlation is completely reversed by the tt pulses (see 
the discussion in Section 0, -/Vfljp.A (i) is also reversed 
and returns to zero at each spin echo time. Therefore, 
^flip.A {t) does not accumulate in the pulse controlled 
dynamicSi^S and we just need to look at the maximum 
value of -/Vflip,A {t) between echoes. 

For single system FID and ensemble spin echo calcula- 
tion presented in this paper, the sufficient condition for 
PC A to be valid is N'l^^ p^{Th) < by noticing that 

T2 « min(T2,A,T2,B) and Th = y/2T2,B- From Eq. (O, 
we have, 

Nm^MTn) - N-^'/^^l-H-^A', (D3) 

Therefore, we obtain a lower bound on the dot size N 
for the vahdity of PCA: A^^ ^ A^n~%-'^. For GaAs 
fluctuation dot in a field of 10 T, the above condition is 
well satisfied for A^ > 10''. 

In contrast to the non-local pair dynamics, the local 
pair dynamics is not reversed under the influence of the 



electron spin flip and A^flip,B(i) accumulates all through 
the time. Nonetheless, it turns out that for all scenar- 
ios of interest including the evolution under the pulse 
sequence control as in22i, we have, 

A^flip.B {t) < Bit' ^ Afsbh' ^ Nbh' (D4) 
fc 

where Afs ~ A' is the number of local nuclear spin pairs. 
For the relevance of the single system FID and ensemble 
spin echo calculation, we shall examine 

Afflip.BCTff) « N^^HA-^ (D5) 

and therefore, the condition (Affiip,B)^ <C A^ sets an upper 
bound on the quantimr dot size A^: N'b'A~' <C 1. In 
GaAs, the above condition is well satisfied for quantum 
dot size Af < 10^ 

To conclude, for nuclear pair-correlations to dominate 
over higher order correlations in the relevant timescale, 
local nuclear pair-dynamics driven by the intrinsic cou- 
plings imposes an upper bound on A^ while non-local 
nuclear pair-dynamics driven by the extrinsic coupling 
imposes a lower bound. This mesoscopic regime covers 
quantum dots of all practical size. Within this meso- 
scopic regime, higher order nuclear correlations are well 
negligible under the timescale of interest. The error es- 
timation is based on characterizing the difference in the 
Hilbert space structure of the exact dynamics and that of 
the pseudo-spin model and assuming this difference has a 
full influence on the electron spin coherence calculation. 
Therefore, the bound of Eq. IClip is not necessarily a 
tight bound. 
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